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Abstract 

All current models of our Universe based on General Relativity have in common 
that space is presently in a state of expansion. In this expository paper we ad- 
dress the question of whether, and to what extent, this expansion influences the 
dynamics on small scales (as compared to cosmological ones), particularly in our 
Solar System. Here our reference order of magnitude for any effect is given by 
the apparent anomalous acceleration of the Pioneer 10 and 1 1 spacecrafts, which 
is of the order of lO^^m/s'^. We distinguish between dynamical and kinematical 
effects and critically review the status of both as presented in the current Utera- 
ture. We argue that in the Solar System dynamical effects can be safely modeled 
by suitably improved Newtonian equations, showing that effects do exist but are 
smaller by many orders of magnitude compared to our Pioneer reference. On 
the other hand, the kinematical effects need a proper relativistic treatment and 
have been argued by others to give rise to an additional acceleration of the order 
He, where H is the Hubble parameter and c is the velocity of light. This simple 
and suggestive expression is intriguingly close to the anomalous Pioneer accel- 
eration. We reanalyzed this argument and found a discrepancy by a factor of 
{v/cY, which strongly suppresses the alleged i/c-effect for the Pioneer space- 
crafts by 13 orders of magnitude. We conclude with a general discussion which 
stresses the fundamental importance to understand precisely, i.e. within the full 
dynamical theory (General Relativity), the back-reaction effects of local inho- 
mogeneities for our interpretation of cosmological data, a task which is not yet 
fully accomplished. Finally, a structured literature list of more than 80 references 
gives an overview over the relevant publications (known to us). 



1 



Contents 



1 Introduction 

2 The Newtonian approach 

2. 1 The two-body problem in an expanding universe 

2.2 Specifying the initial-value problem 

2.3 Discussion of the reduced effective potential 

3 Fully relativistic treatment for gravitationally-bounded systems: 
The McVittie model 

3. 1 Interpretation of the McVittie metric 

3.2 Motion of a test particle in the McVittie spacetime 

4 Fully relativistic treatment for electromagnetically-bounded systems 
and the argument of Dicke and Peebles 

4. 1 The argument of Dicke and Peebles 

4.2 Equations of motion in the physical coordinates 

5 Kinematical effects 

5. 1 Local Einstein-simultaneity in general spacetimes 

5.2 Application to isotropic cosmological metrics 

6 Summary and outlook 

A Additional material 

A.l A simple estimate of the dynamical effect of cosmological expansion . 

A.2 The Schiicking radius for different astronomical scales 

A. 3 Some astronomical and cosmological data 



1 Introduction 



The overall theme of this paper concerns the question of whether the global cosmolog- 
ical expansion has any influence on the local dynamics and kinematics within the Solar 
System. Despite many efforts in the past, this problem is still debated controversially 
in the current hterature and hardly anything is written about it in standard textbooks. 
This is rather strange as the question seems to be of obvious interest. Hence there is 
room for speculations that such an influence might exist and be detectable with current 
experimental means. It has even be suggested that it be (partly) responsible for the 
apparently anomalous acceleration of the Pioneer spacecrafts, the so-called 'Pioneer- 
Anomaly' [74, 75, 78, 79, 80, 81], henceforth abbreviated by PA. 

Existing investigations in this direction arrive at partially conflicting conclusions. To 
resolve this issues at a fundamental level one may think to start from an investigation 
of the embedding-problem for the Solar System. From the viewpoint of General Rela- 
tivity this means to at least find a solution to Einstein's field equations that match the 
gravitational field (i.e. the spacetime metric g^^) of the Sun with that of the inner- 
galactic neighbourhood. In a first approximation the region exterior to some radius 
outside the Sun may be modelled by a spherically symmetric (around the Sun's center) 
constant dust-Uke mass distribution (i.e. matching an exterior Schwarzschild solution 
to a Friedman-Robertson- Walker solution). Further refinements may then take into ac- 
count the structured nature of the cosmological mass distribution as well as a matching 
of the gravitational field of the Galaxy to that of its cosmological environment. 

Analytical solutions are known for some embedding problems under various special 
assumptions concerning symmetries and the matter content of the cosmological envi- 
ronment: 

A. Spherical symmetry, matching Schwarzschild to a Friedmann-Robertson- 
Walker (FRW) universe with pressureless matter and zero cosmological constant 
A. [7, 34]. 

B. Obvious generaUzations of A to many, non-overlapping regions — so-called 
'swiss-cheese models'. 

C. Generalizations of A,B for A 7^ [1]. 

D. Generahzations of A to spherically symmetric but inhomogeneous Lemaitre- 
Tolman-Bondi cosmological backgrounds [3]. 

In all these approaches, which are based on the original idea of Einstein & Straus [7], 
the matching is between a strict vacuum solution (Schwarzschild) and some cosmo- 
logical model. Hence, within the matching radius, the solution is strictly of the 
Schwarzschild type, so that there will be absolutely no dynamical influence of cos- 
mological expansion on local dynamical processes by construction. The only relevant 
quantity to be determined is the matching radius, rs, (sometimes called the 'Schiicking 
radius') as a function of the central mass and the time. This can be done quite straight- 
forwardly and it turns out to be expressible in the following form: 

Vol(r5) ■p = M. (1) 
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Here Vol(r5) is the 3 -dimensional volume within a sphere of radius rs in the cos- 
mological background geometry and p is the (spatially constant but time dependent) 
cosmological mass-density. For example, for flat or nearly flat geometries we have 
Vo^r^) = 47rr|/3. Thus the defining equation for the matching radius has the fol- 
lowing simple interpretation: if you want to embed a point mass M into a cosmolog- 
ical background of mass density p you have to place it at the centre of an excised ball 
whose mass (represented by the left hand side of Q) is just M. This is just the obvious 
dynamical condition that the gravitational pull the central mass M exerts on the am- 
bient cosmological masses is just the same as that of the original homogeneous mass 
distribution within the ball. The deeper reason for why this argument works is that the 
external gravitational field of a spherically symmetric mass distribution is independent 
of its radial distribution, i.e. just depending on the total mass. This well known fact 
from Newton's theory remains true in General Relativity, as one readily sees by recall- 
ing the uniqueness of the exterior Schwarzschild solution. We may hence simply think 
of the cosmological mass inside a sphere of radius rs as being squashed to a point (or 
black hole) at its centre, without affecting the dynamics outside the radius rs- Note 
that for expanding universes the radius rs expands with it, that is, it is co-moving with 
the cosmological matter. 

At first sight this simple solution may appear as a convincing argument against any 
influence of cosmic expansion on the scale of our Solar System, as rs certainly extends 
far beyond it; too far in fact ! If M is the mass of the Sun then rs turns out to be about 
175 pc, which is more than a factor of 100 larger than the distance to our next star and 
more than factor of 50 larger than the average distance of stars in our Galaxy. But this 
means that the Einstein-Straus solution is totally inappropriate as a model for the Solar 
System's neighbourhood. This changes as one goes to larger scales, beyond that of 
galaxy clusters; see Section IaT2| 

The Einstein-Straus solution may also be criticized on theoretical grounds, an obvious 
one being its dynamical instability. Slight perturbations of the matching radius to 
larger radii will let it increase without bound, slight perturbations to smaller radii will 
let it collapse. This can be proven formally (e.g. [22], Ch. 3) but it is also rather 
obvious, since rs is defined by the equal and opposite gravitational pull of the central 
mass on one side and the cosmological mass on the other. Both pulls increase as 
one moves towards their side, so that the equilibrium position must correspond to a 
local maximum of the gravitational potential. Another criticism of the Einstein-Straus 
solution concerns the severe restrictions under which it may be generalized to non 
spherically-symmetric situations; see e.g. [35, 26, 27, 28]. 

We conclude from all this that we cannot expect much useful insight, as regards practi- 
cally relevant dynamical effects within the Solar System, from further studies of mod- 
els based on the Einstein-Straus matching idea. ^ Rather, we shall proceed in the fol- 
lowing steps: 

I. Discuss an improved Newtonian model including a cosmological expansion 
term. This we did and the results are given in Section |2l below. Our discussion 
complements [48] which just makes a perturbative analysis, thereby missing all 

' However, as stressed in Section|6| the matching problem certainly is important on larger scales. 
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orbits which are unstable under cosmological expansion (which do exist). In 
this respect it follows a very similar strategy as proposed in the recent paper by 
Price [65] (the basic idea of which goes back at least to Pachner's work [63, 64]), 
though we think that there are also useful differences. We also supply quantita- 
tive estimates and clarify that the improved Newtonian equations of motion aie 
written in terms of the right coordinates (non-rotating and metrically normal- 
ized). The purpose of this model is to develop a good physical intuition for the 
qualitative as well as quantitative features of any dynamical effects involved. 

II. Eventually the Newtonian model just mentioned has to be understood as a lim- 
iting case of a genuinely relativistic treatment. For the gravitational case this is 
done in Section|3l where we employ the McVittie metric to model a spherically 
symmetric mass embedded in a A; = FRW universe. The geodesic equation 
is then, in a suitable limit, shown to lead to the improved Newtonian model 
discussed above. The same holds for the electromagnetic case, as we show in 
Section 131 There we take a slight detour to also reconsider a classic argument 
by Dicke & Peebles [49], which allegedly shows the absence of any relevant 
dynamical effect of global expansion. Its original form only involved the dy- 
namical action principle together with some simple scaling argument. Since this 
reference is one of the most frequently cited in this field, and since the simplicity 
of the argument (which hardly involves any real analysis) is definitely deceptive, 
we give an independent treatment that makes no use of any scaling behaviour of 
physical quantities other than spatial lengths and times. Our treatment also re- 
veals that the original argument by Dicke & Peebles is insufficient to discuss 
leading order effects of cosmological expansion. It is therefore also ineffective 
in its attempt to contradict Pachner [63, 64]. 

III. Neither the improved Newtonian model nor other general dynamical arguments 
make any statement about possible kinematical effects, i.e. effects in connection 
with measurements of spatial distances and time durations in a cosmological en- 
vironment whose geometry changes with time. This is an important issue since 
tracking a spacecraft means to map out its 'trajectory', which basically means 
to determine its simultaneous spatial distance to the observer at given observer 
times. But we know from General Relativity that the concepts of 'simultane- 
ity' and 'spatial distance' are not uniquely defined. This fact needs to be taken 
due care of when analytical expressions for trajectories, e.g. solutions to the 
equations of motion in some arbitrarily chosen coordinate system, are compared 
with experimental findings. In those situations it is likely that different kinemat- 
ical notions of simultaneity and distance are involved which need to be properly 
transformed into each other before being compared. For example, these trans- 
formations can result in additional acceleration terms which have been claimed 
in the literature to be directly relevant to the PA; see [70, 69, 73, 72, 71, 68]. 
We will confirm the existence of such effects in principle, but are in essential 
disagreement concerning their relevance in practice. We think that they have 
been overestimated by about 13 orders of magnitude. The details will be given 
in Section 121 

IV. Finally we made a systematic scan of the literature on the subject. The papers 
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found to be relevant are listed in the bibliography at the end, which we subdi- 
vided into four sections: . 

1. Papers deaUng with the proper matching problem in General Relativity. 

2. Papers deahng generally with the influence of the global cosmological ex- 
pansion on local dynamics, irrespectively of whether they work within an 
improved Newtonian setting or in full General Relativity. 

3. Papers discussing tentative explanations of the PA by means of gravity, 
mostly by referring to kinematical effects of space-time measurements in 
time dependent background geometries. 

4. Measurements of the PA. These were just for our own instruction and are 
Usted for completeness. 

We beUeve one can give a fair estimation on the irrelevance of the dynamical effects 
in question. Given the weakness of the gravitational fields involved, estimations by 
Newtonian methods should give reliable figures of orders of magnitude for the motion 
of ordinary matter. Kinematical effects based on the equation for light propagation in 
an expanding background also turn out to be negligible, in contrast to some claims in 
the literature, which we think can be straightened out. We end by suggesting possible 
routes for further research. 

2 The Newtonian approach 

In order to gain intuition we consider a simple bounded system, say an atom or a 
planetary system, immersed in an expanding cosmos. We ask for the effects of this 
expansion on our local system. Does our system expand with the cosmos? Does it 
expand only partially? Or does it not expand at all? 

2.1 The two-body problem in an expanding universe 

Take a two-body problem with a attractive force between them. For simplicity we 
think of one mass as being much smaller than the other one (this is inessential). This 
can e.g. be a system consisting of two galaxies, a star and a planet (or spacecraft), or a 
(classical) atom given by an electron orbiting around a proton. We think this system as 
being immersed in an expanding universe and we model the effect of the cosmological 
expansion by adding to the attraction term an extra term coming from the Hubble law 
r = Hr. Here H := a/a is the Hubble parameter, a{t) the cosmological scale factor, 
and r the distance - as measured in the surface of constant cosmological time t - of 
two objects that follow the Hubble flow (cosmological expansion). The acceleration 
that results from the Hubble law is 

cosm.acc. ~ • (2) 

a 

Note that, in the sense of General Relativity, a body that is co-moving with the cosmo- 
logical expansion is moving on an inertial trajectory, i.e. it moves force free. Forces in 
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the Newtonian sense are now the cause for deviations form the co-moving acceleration 
described by This suggests that in Newton's law, mf = F, we have to replace r 
by r — {a/a)r. This can be justified rigorously by using the equation of geodesic devi- 
ation in General Relativity. In order to do this one must make sure that the Newtonian 
equations of motion are written in appropriate coordinates. That is, they must refer 
to a (locally) non-rotating frame and directly give the spatial geodesic distance. This 
is achieved by using Fermi normal coordinates along the worldline of a geodesically 
moving observer — in our case e.g. the Sun or the proton — , as correctly emphasized in 
[48]. The equation of geodesic deviation in these coordinates now gives the variation 
of the spatial geodesic distance to a neighbouring geodesically moving object — in our 
case e.g. the planet (or spacecraft) or electron. It reads^ 

-^ + R\i^x' = Q. (3) 

Here the x'^ are the spatial non-rotating normal coordinates whose values directly refer 
to the proper spatial distance. In these coordinates we further have [48] 

R\io = -8la/a (4) 

on the worldline of the first observer, where the overdot refers to differentiation with 
respect to the cosmological time, which reduces to the eigentime along the observer's 
worldline. 

Neglecting large velocity effects (i.e. terms quadratic or higher order inv/c) we can 
now write down the equation of motion for the familiar two-body problem. After 
specification of a scale function a{t), we get two ODEs for the variables (r, Lp), which 
describe the position^ of the orbiting body with respect to the central one: 

.. L'^ C a 

r = ^ ^ + -r (5a) 

r'^ip = L. (5b) 

These are the (a/a)-improved Newtonian equations of motion for the two-body prob- 
lem, where L represents the (conserved) angular momentum of the planet (or elec- 
tron) per unit mass and C the strength of the attractive force. In the gravitational case 
C = GM, where M is the mass of the central body, and in the electromagnetic case 
C = \Qq\/47reom (Sl-unit), where Qq is the product of the two charges, m is the 
electron mass, and eo the vacuum permittivity. In Sections|3]and|4]we will show how 
to obtain ^ in appropriate limits from the full general relativistic treatments. 

We now wish to study the effect the a term has on the unperturbed Kepler orbits. We 
first make the obvious remark that this term results from the acceleration and not just 
the expansion of the universe. 

^ By construction of the coordinates, tiie Qiristoffel symbols F^^ vanish along the worldhne of the 
first observer. Since this worldline is geodesic, Fermi-Walker transportation just reduces to parallel 
transportation. This gives a non-rotating reference frame that can be physically realized by gyros 
taken along the worldline. 

^ Recall that 'position' refers to Fermi normal coordinates, i.e. r is the radial geodesic distance to the 
observer at r = 0. 
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We first remark that, in the concrete physical cases of interest, the time dependence of 
this term is neghgible to a very good approximation. Indeed, putting / := a/ a, the 
relative time variation of the coefficient of r in (IJli is ///. For an exponential scale 
function a{t) oc exp(Ai) (A-dominated universe) this vanishes, and for a power law 
a{t) oc (for example matter-, or radiation-dominated universes) this is —2H/\, and 
hence of the order of the inverse age of the universe. If we consider a planet in the 
Solar System, the relevant time scale of the problem is the period of its orbit around 
the Sun. The relative error in the disturbance, when treating the factor a/a as constant 
during an orbit, is hence smaller than 10^^ for the planets in the Solar System. For 
atoms it is much smaller, of course. Henceforth we shall neglect this time-dependence 
of©. 

Keeping this in mind we set from now on a/a = const =: A. Taking the actual value 
one can write A = —qoHq, where the index zero means 'today' and q stands for the 
cosmological deceleration parameter (defined by q := —a/{H'^a)). Since the force is 
time independent, we can immediately integrate (Bat and get 

lr2 + y(r)=i?, (6) 

where the effective potential is 



2.2 Specifying the initial-value problem 

For ^ and dSbl we have to specify initial conditions (r, r, y?, ip){tQ) = (ro, vq, c/Jq, wq) 
at the initial time t^. To study the solutions of the above equation for r one has to 
look at the effective potential. For this purpose it is very convenient to introduce a 
length scale and a time scale that naturally arise in the problem. The length scale is 
defined as the radius at which the acceleration due to the cosmological expansion has 
the same magnitude as the gravitational (or electromagnetical) attraction. This happens 
precisely at the critical radius 

/ (7 \ 1/3 

[\A\) ■ 

For r < r^rit the gravitational (or electromagnetical) attraction dominates, whereas for 
r > TcHt the effect of the cosmological expansion is the dominant one. 



Intermezzo: Expressing the critical radius in terms of cosmological parameters 

We briefly wish to point out how to express the critical radius in terms of the cosmo- 
logical parameters. For this we write: 

GM f M^/^ 

120pc, (9) 
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where we have used the current values qq = —1/2 and /iq = 0.7. It is interesting to 
note that in the case of zero cosmological constant and pressureless matter we recover 
the Schiicking gluing radius [34] (compare ©): 



rs = (tt^- ) ■ (10) 



1/3 

V(4/3)7r/>," 

To see this, just recall that for a pressureless matter we have go = (l/2)J^m — ^^A- 
Then, for a vanishing cosmological constant, and using the definition il^ = Pm • 
8ttG/ (S/Zq ), one gets the above equations immediately. 

In the electromagnetic case, for a proton-electron system, 

^'^'-30AU, (11) 



^47reom|go|-f^o 
which is about as big as the Neptune orbit! 



Tk -.= 271 {^] . (12) 



Back to the initial-value problem 

The time scale we define is the period with respect to the unperturbed Kepler orbit (a 
solution to the above problem for ^4 = 0) of semi-major axis tq. By Kepler's third law 
it is given by 

,3\ 1/2 
0. 

It is convenient to introduce two dimensionless parameters which essentially encode 
the initial conditions (ro, wq)- 

a := sign(^) ( ^) = A^^ . (14) 

For close to Keplerian orbits A is close to one. For reasonably sized orbits a is close 
to zero. For example, in the Solar System, where rp < 100 AU, one has |q| < 10~^^. 
For an atom whose radius is smaller than 10^ Bohr-radii we have \a\ < 10~^^. 

Defining 

x{t) := r{t)/ro , (15) 
equations ^ and dSbt can now be written as 

h^ + {2TT/TK)^vx,a{x)=e (16) 
x'^ip = loq , (17) 

where e := E/r^ now plays the role of the energy-constant and where the reduced 
2-parameter effective potential vx^a is given by 

A 1 a 2 
vx.a{x) ■= TT^ t:x ■ (18) 
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The initial conditions now read 

(x, X, V?, ip){to) = (1, vo/ro, 9?o, Wo) . (19) 

The point of introducing the dimensionless variables is that the three initial par ameters 
(L, C, A) of the effective potential could be reduced to two: A and a). This will be 
convenient in the discussion of the potential. 

2.3 Discussion of the reduced effective potential 




Figure 1: The figure shows the effective potential fA,«, for circular orbits, where A = 
1 — a, for some values of a. The initial conditions are x = 1 and x = (see dTSl ). 
At X = 1 the potential has an extremum, which for a < 1/4 is a local minimum 
corresponding to stable circular orbits. For 1/4 < a < 1 these become unstable. 

Circular orbits correspond to extrema of the effective potential 0. Expressed in terms 
of the dimensionless variables this is equivalent to v'^^{l) = — A + 1 — a = 0. 
By its very definition (fT3b . A is always nonnegative, implying a < 1. For negative a 
(decelerating case) this is always satisfied. On the contrary, for positive a (accelerating 
case), this implies the existence of a critical radius 

ro < r.rit (20) 

beyond which no circular orbit exists. 
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These orbits are stable if the considered extremum is a true minimum, i.e. if the second 
derivative of the potential evaluated at the critical value is positive. Now, v"^{l) = 
3A — 2 — a = l — 4a, showing stability for a < 1/4 and instability for a > 1/4. 

Expressing this in physical quantities we can summarize the situation as follows: in the 
decelerating case (i.e. for negative a or, equivalently, for negative A) stable circular 
orbits exist for every radius tq; one just has to increase the angular velocity according 
to (1211) . On the contrary, in the accelerating case (i.e. for positive a, or, equivalently, 
for positive A), we have three regions: 

• ^0 < '"ub := (l/4)^/'^reHt ~ 0.63r„it, where circular orbits exist and are stable.^ 

• ^^.b ^ fQ < r<,Ht, where circular orbits exist but are unstable. 

• ro > r^rit, where no circular orbits exist. 



Generally, there exist no bounded orbits that extend beyond the critical radius r^rit, the 
reason being simply that there is no r > r„it where V'{r) > 0. Bigger systems will 
just be slowly pulled apart by the cosmological acceleration and approximately move 
with the Hubble flow at later tunes.^ 

Turning back to the case of circular orbits, we now express the condition for an extrema 
derived above, A = 1 — a, in terms of the physical quantities, which leads to 



ojo = (27r/Tx)Vl-sign(yl)(ro/r„,,)3 . (21) 

This equation says that, in order to get a circular orbit, our planet, or electron, must 
have a smaller or bigger angular velocity according to the universe expanding in an 
accelerating or decelerating fashion respectively. This is just what one would expect, 
since the effect of a cosmological 'pulling apart' or 'pushing together' must be com- 
pensated by a smaller or larger centrifugal forces respectively, as compared to the 
Keplerian case. Equation dlTl represents a modification of the third Kepler law due 
to the cosmological expansion. In principle this is measurable, but it is an effect of 
order {ro/r^^n)^ and hence very small indeed; e.g. smaller than 10^^^ for a planet in 
the Solar System. 

Instead of adjusting the initial angular velocity as in d2n i. we can ask how one has 
to modify tq in order to get a circular orbit with the angular velocity loq = 2it/Tk- 
This is equivalent to searching the minimum of the effective potential (I18t for A = 1. 
This condition leads to the fourth order equation ax'^ — x + 1 = with respect to 
X. Its solutions can be exactly written down using Ferrari's formula, though this is 
not illuminating. For our purposes it is more convenient to solve it approximatively, 
treating a as a small perturbation. Inserting the ansatz Xmin = cq + cia + O(a^) we 
get Co = ci = 1. This is really a minimum since v'l^{xmin) = 1 + 0{a) > 0. Hence 
we have 

rn.in = ro ( 1 + sign(A) ( —\ + oUr^/r^^.^f ) ) (22) 



crit 



'ub' stands for 'upper bound for stable circular orbits' 
^ This genuine non-perturbative behaviour was not seen in the perturbation analysis performed in [48]. 
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This tells us that in the accelerating (decelerating) case the radii of the circular orbits 
with loq = 2it/Tk becomes bigger (smaller), again according to expectation. As an 
example, the deviation in the radius for an hypothetical spacecraft orbiting around the 
Sun at 100 AU would be just of the order of 1 mm. Since it grows with the fourth 
power of the distance, the deviation at 1000 AU would be be of the order of 10 meters. 

3 Fully relativistic treatment for gravitationally-bounded 
systems: The McVittie model 

In physics we are hardly ever in the position to mathematically rigorously model phys- 
ically realistic scenarios. Usually we are at best either able to provide approximate 
solutions for realistic models or exact solutions for approximate models, and in most 
cases approximations are made on both sides. The art of physics then precisely con- 
sists in finding the right mixture in each given case. However, in this process our 
intuition usually strongly rests on the existence of at least some 'nearby' exact solu- 
tions. Accordingly, in this section we seek to find exact solutions in General Relativity 
that, with some degree of physical approximation, model a spherically symmetric body 
immersed in an expanding universe. 

There are basically two ways to proceed, which could be described as the 'gluing'- and 
the 'melting' -way respectively. In the first and simpler approach one constructs a new 
solution to Einstein's equations by suitably gluing together two known solutions, one 
corresponding to the star (i.e. Schwarzschild, if the star has negligible angular momen- 
tum), the other to a homogeneous universe (i.e. FRW). The resulting spacetime is then 
divided into two distinct regions, whose interiors are locally isometric to the original 
solutions. The Einstein-Straus-Schiicking vacuole [7, 34] clearly belongs to this class, 
as well as its generalizations [1, 3]. The advantage of this gluing-approach is its rela- 
tive analytic simplicity, since the solutions to be matched are already known. Einstein's 
equations merely reduce to the junction conditions along their common seam. Its dis- 
advantage is that this gluing only works under those very special conditions which 
allow the glued solutions to locally persist exactly, and these conditions are likely to 
be physically unrealistic. 

In the second approach one considers genuinely new solutions of Einstein's equa- 
tions which only approximately resemble the spacetimes of an isolated star or a ho- 
mogeneous universe at small and large spatial distances respectively. This 'melting- 
together' is the more flexible approach which therefore allows to model physically 
more realistic situations. Needless to say that it also tends to be analytically more 
complicated and that a physical interpretation is often not at all obvious. The solutions 
of McVittie [59] and also that of Gautreau [14] fall under this class. 

Our goal is to accurately model the Solar System in the currently expanding spatially- 
flat universe. We already argued in Section ^that the Einstein-Straus-Schiicking vac- 
uole model unfortunately does not apply to the Solar System since the matching radius 
would be much too big (see also Appendix IA.2I) . The model of Gautreau [14] is much 
harder to judge. Its physical and mathematical assumptions are rather implicit and not 
easy to interpret as regards their suitability for the problem at hand. For example, it 
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assumes the cosmic matter to move geodesically outside the central body, but at the 
same time also assumes an equation of state in the form p{p). For p 7^ this seems 
contradictory in a genuinely non-homogeneous situation since pressure gradients will 
necessarily result in deviations from geodesic motions (see [38]). Assuming p = 
(which implies a motion of cosmic matter with non-vanishing shear) Gautreau finds 
in [14] that orbits will spiral into the central mass simply because there is a net influx 
of cosmic matter and hence an attracting source of increasing strength. This is not 
really the kind of effect we are interested in here. 

Among the models discussed in the literature the one that is best understood as regards 
its analytical structure as well as its physical assumptions is that of McVittie [59] for 

= 0. It is therefore, in our opinion, the natural candidate to consider first when 
modeling systems like the Solar System in our expanding Universe. As emphasized 
above, this is not to say that this model is to be considered realistic in all its detailed 
aspects, but at least we do have some fairly good control over the assumptions it is 
based on thanks to the carefully analysis by Nolan [30, 31, 32]. For example, there is a 
somewhat unrealistic behaviour of the McVittie spacetime and the matter in it near the 
singularity at r = m/2, as discussed below. But at larger radial distances the 'flat' (i.e. 

= 0) McVittie model may well give a useful description of the exterior region of a 
central object in an expanding spatially-flat universe, at least in the region where the 
radius is much larger (in geometric units) than the central mass (to be defined below). 
For planetary motion and spacecraft navigation in the Solar System this is certainly 
the case, since the ratio of the central mass to the orbital radius is of the order of 
1.5km/lAU = 10"®. 

In this section we briefly review the McVittie spacetime and look at the geodesic equa- 
tion in it, showing that it reduces to ^ in an appropriate weak-field and slow-motion 
hmit. This provides another more solid justification for the Newtonian approach we 
carried out in Section|2l 

The flat McVittie model (from now on to be simply referred to as the McVittie model) 
is characterized by two inputs: First, one makes the following ansatz for the metric 
which represents an obvious attempt to melt together the Schwarzschild metric (in 
spatially isotropic coordinates) with the spatially flat FRW metric (I37t : 

\l + m{t)/2r J \ 2r J ^ ' ^ ' 

Here dil^ = dO"^ + sin^ dip^ and the two time-dependent functions m and a are 
to be determined. It is rotationally symmetric with the spheres of constant radius 
being the orbits of the rotation group. Second, it is assumed that the ideal fluid (with 
isotropic pressure) representing the cosmological matter moves along integral curves 
of the vector field d/dt. Note that this vector field is not geodesic (unlike in the 
Gautreau model). Moreover, this model contains the implicit assumption that the fluid 
motion is shearless (cf. Chapter 16 of [36]). 

Einstein's equations together with the equation of state determine the four functions 
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m{t),a{t), p{t, r), and p{t, r). The former are equivalent to 



[am) = 



(24a) 




2 



Svr/) 



(24b) 



Svrp 



) 



( 



1 + m/2r 



1 — m/2r 



) 



(24c) 



Here we already used the first one in order to express the derivatives of m in terms of 
a and its derivatives. The first equation can be immediately integrated: 



where mo is an integration constant. Below we will show that this integration constant 
is to be interpreted as the mass of the central particle. We will call the metric d23b 
together with condition d25t the McVittie metric. 

We note that the equation of state must be necessarily space dependent. This follows 
directly from the equations (I24bl) and (124c t . which imply that the density depends 
on the time coordinate only whereas the pressure depends on both time and space 
coordinates. Formally the system S24i can be looked upon in two ways: either one 
prescribes an equation of state and deduces from (I24bl) . (I24cl) . and d25t a second-order 
differential equation for the scale factor a, or one specifies a{t) and deduces the matter 
density, the pressure, and hence the equation of state. 

As special cases of d24t we remark that if either a or m are time independent (I24at 
implies that both must be time independent. This, in turn, implies that density and pres- 
sure vanish everywhere, resulting in the Schwarzschild solution in spatially isotropic 
coordinates. If we choose the equation of state to be that of pressureless dust, i.e. 
p = 0, we get either the Schwarzschild solution or the dust-filled FRW universe. This 
follows from (I24cl) . where we must distinguish between two cases: d/a can only be 
constant if it is zero, hence resulting in the Schwarzschild solution. If d/a is not con- 
stant, equation (I24ct (with d25t ) implies, after a partial differentiation with respect to r, 
that niQ = 0. This gives the homogeneous and isotropic dust-filled FRW universe. An- 
other possible equation of state is that of a cosmological-constant. This choice implies 
constancy of a/ a and hence that the second term on the left hand side of (124c t van- 
ishes. In this way one recovers the Schwarzschild-de Sitter metric in spatially isotropic 
coordinates. 

Finally we also mention some critical aspects of the McVittie model. In fact, unless 
d/a is constant, it has a singularity at r = m/2 where the pressure as well as some 
curvature invariants diverge. The former can be immediately seen from (I24cl) . This 
is clearly a result of the assumption that the fluid moves along the integral curves of 
d/dt, which become lightlike in the limit as r tends to m/2. Their acceleration is 
given by the gradient of the pressure, which diverges in that limit. For a study of the 
singularity at r = m/2 see [31, 32]. 

A related question concerns the global behaviour of the McVittie metric. Each hy- 
persurface of constant time t is a complete Riemannian manifold which besides the 
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rotational symmetry admits a discrete isometry, given in (r, 9, cp) coordinates by 

(Pir, 9, ^) = ([mo/2a(t)]2 r'^ 9 , y,) . (26) 

It corresponds to a reflection at the 2-sphere r = (mo/2a(t)) and shows that the hy- 
persurfaces of constant t can be thought of as two isometric asymptotically-flat pieces 
joined together at the totally geodesic (being a fixed-point set of an isometry) 2-sphere 
r = mo/2a(t), which is minimal. Except for the time-dependent factor a{t), this is 
just like for the slices of constant t in the Schwarzschild metric (the difference being 
that d26l does not extend to an isometry of the spacetime metric unless a = 0). This 
means that the McVittie metric cannot literally be interpreted as corresponding to a 
point particle sitting at r = (r = is in infinite metric distance) in a flat FRW uni- 
verse, just like the Schwarzschild metric does not correspond to a point particle sitting 
at r = in Minkowski space. Unfortunately, McVittie seems to have interpreted his 
solution in this fashion [59] which even until recently gave rise to some confusion in 
the literature (e.g. [14, 37, 12]). A clarification was given by Nolan [31]. Another 
important issue is whether the cosmological matter satisfies some energy condition. 
For a discussion about this topic we refer to [30, 31]. 

3.1 Interpretation of the McVittie metric 

We shall now present some arguments which justify calling McVittie's metric a model 
for a localized mass immersed in a flat FRW background. Here we basically fol- 
low [30]. As is well known, it is generally not possible in General Relativity to assign 
a definite mass (or energy) to a local bounded region of space (quasi-local mass). Phys- 
ically sensible definitions of such a concept of quasi-local mass exist only in favourable 
and special circumstances, one of them being spherical symmetry. In this case the so- 
called Misner-Sharp energy is often employed (e.g. [29, 5, 6, 4, 17, 18]). It allows 
to assign an energy content to the interior region of any two-sphere of symmetry (i.e. 
an orbit of 50(3)). For the McVittie metric ((|23ll with (|25l)) the Misner-Sharp energy 
takes the simple and intuitively appealing form: 

E^sig; R, t) = ^irR^pit) + mo (27) 

where henceforth we denote by R the 'areal radius' defined by ilSl . This shows that 
the total energy is given by the sum of the cosmological matter contribution and the 
central mass, where the mass of the central object is given by niQ. 

Another useful definition of quasi-local mass is that of Hawking [16]. According to 
this definition the energy contained in the region enclosed by a spatial two-sphere S 
is given by a surface integral over S, whose integrand is essentially the sum of certain 
distinguished components of the Ricci and Weyl tensors, representing the contributions 
of matter and the gravitational field respectively. Applied to the McVittie metric the 
latter takes the value ttiq for any 2-sphere outside of and enclosing R = 2mQ. 
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3.2 Motion of a test particle in the McVittie spacetime 



We are interested in the motion of a test particle (idealizing a planet or a spacecraft) 
in McVittie's spacetime. In [59] McVittie concluded within a slow-motion and weak- 
field approximation that Keplerian orbits do not expand as measured with the 'cos- 
mological geodesic radius' = a{t)r. Later Pachner [63] and Noerdlinger & Pet- 
rosian [62] argued for the presence of the acceleration term ^ proportional to d/a 
within this approximation scheme, hence arriving at d5at . In the following we shall 
show how to arrive at dSal i from the exact geodesic equation of the McVittie metric by 
making clear the approximations involved. In order to compare our calculation with 
similar ones in the recent literature (i.e. [43, 41])^ we will work with the so-called 
'areal radius'. It corresponds to a function that can be geometrically characterized 
on any spherically symmetric spacetime by taking the square root of the area of the 
50(3)-orbit through the considered point divided by 47r. Hence it is the same as the 
square root of the modulus of the coefficient of the angular part of the metric. For the 
McVittie metric this reads 

R{t,r)= + a{t)r. (28) 



2a{t)r 

Note that for fixed t the map r i-^ R{t, r) is 2-to-l and that R > 2mo, where R = 2mo 
corresponds to r = mo/2a. Hence we restrict the coordinate transformation (I28t to 
the region r > mo /2a where it becomes a diffeomorphism onto the region R > 2mo. 

Reintroducing factors of c, McVittie's metric assumes the (non-diagonal) form in the 
region R > 2mo (i.e. r > mo/2a(t)) 

where we put 

f{R) := 1 - ^ , (30) 

H{t):=-{t). (31) 
a 

The region R < 2mo was investigated in [32]. 

The equations for a timelike geodesic (i.e. parameterized with respect to eigentime) 
r 1-^ z^(r) with g{z,z) = follows via variational principle from the Lagrangian 
C{z,z) = {l/2)g^iy{z)z'^z'^ . Spherical symmetry implies conservation of angular 
momentum. Hence we may choose the particle orbit to lie in the equatorial plane 
9 = 7r/2. The constant modulus of angular momentum is 

i?V = L . (32) 

The remaining two equations are then coupled second-order ODEs for f(r) and R{t). 
However, we may replace the first one by its first integral that results from g{z, z) = (?: 

^ The paper [43] contains a derivation of the effect of cosmological expansion on the periastron preces- 
sion and eccentricity change in the case where the Hubble parameter H := a/ a is constant. 
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2^ 

V c J r " ' ./fiR) f{R) 



m) - (^) ] cH^ + ™^cti. - ^ - iL/Rf = . (33) 



The remaining radial equation is given by 

+ ^/(-R) 04b) 
- + mtf[l - ^ - (^)')) (340 



(34d) 



mo f H{t)R\ 
~R 



2^ 



Recall that rn-o = GM/c^, where M is the mass of the central star (the Sun in our 
case) in standard units (kg). 

Equations (I33l34t are exact. We are interested in orbits of slow-motion (compared 
with the speed of light) in the region where 

2mo =:Rs<^R<^Rh ■■= c/H . (35) 

The latter condition clearly covers all situations of practical applicability in the Solar 
System, since the Schwarzschild radius Rs of the Sun is about 3 km = 2 • 10"*^ AU 
and the 'Hubble radius' Rh is about 13.7 • 10^ ly = 8.7 • 10^'' AU. 

The approximation now consists in considering small perturbations of Keplerian orbits. 
Let r be a typical timescale of the problem, like the period for closed orbits or else 
R/v with V a typical velocity. The expansion is then with respect to the following two 
parameters: 

- PS ( — -] (slow-motion and weak-field) , (36a) 
c \ R J 

£2 ~ HT (small ratio of chai^acteristic-time to world-age) , (36b) 



In order to make the expression to be approximated dimensionless we multiply d33b 
by and d34b by T"^ / R. Then we expand the right hand sides in powers of the 
parameters (l36l . using the fact that {HR/c) £i£2- From this and (l32l we obtain ^ 
if we keep only terms to zero-order in e-^ and leading (i.e. quadratic) order in 62, where 
we also re-express R as function of t. Note that in this approximation the areal radius 
R is equal to the spatial geodesic distance on the t = const, hypersurfaces. 
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4 Fully relativistic treatment for electromagnetically- 
bounded systems and the argument of Dicke and Peebles 



In this section we show how to arrive at ^ from a fully relativistic treatment of an 
electromagnetically bounded two-body problem in an expanding (spatially flat) uni- 
verse. This implies solving Maxwell's equations in the cosmological background d37l 
for an electric point charge (the proton) and then integrate the Lorentz equations for the 
motion of a particle (electron) in a bound orbit (cf. [45]). Equation ^ then appears 
in an appropriate slow-motion limit. However, in oder to relate this straightforward 
method to a famous argument of Dicke & Peebles, we shall proceed by taking a slight 
detour which makes use of the conformal properties of Maxwell's equations. 



4.1 The argument of Dicke and Peebles 

In reference [49] Dicke & Peebles presented an apparently very general and elegant 
argument, that purports to show the insignificance of any dynamical effect of cosmo- 
logical expansion on a local system that is either bound by electromagnetic or gravita- 
tional forces which should hold true at any scale. Their argument involves a rescaling 
of spacetime coordinates, (t, x) i— > (At, Ax) and certain assumptions on how other 
physical quantities, most prominently mass, behave under such scaling transforma- 
tions. For example, they assume mass to transform like m X'^m. However, their 
argument is really independent of such assumptions, as we shall show below. We work 
from first principles to clearly display all assumptions made. 

We consider the motion of a charged point particle in an electromagnetic field. The 
whole system, i.e. particle plus electromagnetic field, is placed into a cosmological 
FRW-spacetime with flat {k = 0) spatial geometry. The spacetime metric reads 

g = c^ dt^ - a^{t){dr^ + {d9^ + sin^ 6 dip"^)) . (37) 

We introduce conformal time, tc, via 

/■* df' 

tc = fit) := / ^ , (38) 



Jk a{t') 

by means of which we can write (l37l in a conformally flat form, where rj denotes the 
flat Minkowski metric: 

g = al{tc){c^ dtl - dr'^ - [dO'^ + sin^ 9 d^p'^)} . (39) 

^ V ' 

V 

Here we wrote Oc to indicate that we now expressed the expansion parameter a as 
function of tc rather than t, i.e. 

ac:=ao /"^ . (40) 
For example, if a{t) = af^ (0 < n < 1), then we can choose /c = in d38l and have 

f* dt' t^~^ 
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so that 

t = r\Q = [{l-n)at,]'^^'-''\ (42) 

and therefore 

a,{Q = at^/(i-") , where a := [(1 - n^a] ^^^^""^ . (43) 

The electromagnetic field is characterized by the tensor F^,^, comprising electric and 
magnetic fields: 

= r ^-'^ ] . (44) 



In terms of the electi^omagnetic four-vector potential, = {ip/c, — A), one has 

i^M^ = d^A^ - d^A^ = V^A^ - VyA^ , (45) 

so that, as usual, E = — V0 — A. The expression for the four- vector of the Lorentz- 
force of a particle of charge e moving in the field F^^, is eF^aU^, where u is the 
particle's four velocity. 

The equations of motion for the system Particle -i- EM-Field follow from an action 
which is the sum of the action of the particle, the action for its interaction with the 
electromagnetic field, and the action for the free field, all placed in the background 
d37l . Hence we write: 

S = Sp + Si + Sf, (46) 

where 



Sp = —m(?j dr = —mc j \J g{z' , z') dX , (47a) 

Si = -e j^A^dxi" = -e j A^{z{X))z'f' dX 

= - j d^x Af,{x) j dXe 5^^^ {x -z{X)) z'^ , (47b) 

Sf = d"xy^dit^<7^°5"^F^,F„^ = x/ d^^V'^v'^F^uFaf^. (47c) 

Here A is an arbitrary parameter along the worldline z : X ^ z{X) of the particle, and 
z' the derivative dz/dX. The differential of the eigentime along this worldline is 

dr = ^M7;?)dX = ^g^^{z{X))^%dX . (48) 

It is now important to note that 1) the background metric g does not enter (I47bl) and that 
(I47cl) is conformally invariant (in 4 spacetime dimensions only!). Hence the expansion 
factor, a{tc), does not enter these two expressions. For this reason we could write 
(I47ct in terms of the flat Minkowski metric, though it should be kept in mind that 
the time coordinate is now given by conformal time tc- This is not the time read by 
standard clocks that move with the cosmological observers, which rather show the 
cosmological time t (which is the proper time along the geodesic flow of the observer 
field X = d/dt). 
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The situation is rather different for the action (I47al) of the particle. Its variational 
derivative with respect to z(A) is 



SSp \ ^afi,^ z'^'z'f^ d 

^ — -mc ' ^ 



We now introduce the conformal proper time, t^, via 



(49) 



dT, = (1/c) y/7jiz',z')dX = (l/ca) ^/g{z>,z')d\ . (50) 



We denote differentiation with respect to Tc by an overdot, so that e.g. z' / 's/g{z', z') 
z/ca. Using this to replace z' by zy^g{z', z')/ca and also g by a^r] in ( l49t gives 



^ = {^.«^" + (51) 



where we set 



a =: exp((/>/c2) and - -6'^ + ^rj,^ . (52) 

Recalling that 6Sp = J -g^^^dzf^dX = J j^^^6z^^dTc and using d50l . (ISTT i is equiv- 
alent to 



— maiz 



where from now on we agree to raise and lower indices using the Minkowski metric, 
i.e. 7]^,^ = diag(l, —1, —1, —1) in Minkowski inertial coordinates. 

Writing (I47bt in terms of the conformal proper time and taking the variational deriva- 
tive with respect to z{tc) leads to 5Si / 5z^{tc) = —eF^aZ'^, so that 

6S 



ma{z^ + - e . (54) 



The variational derivative of the action with respect to the vector potential A is 

6S 



5,F'^°(x) - e JdT, S{x - z{Te)) zf^iTc) . (55) 



Equations (l54t and (l55l show that the fully dynamical problem can be treated as if it 
were situated in static flat space. The field equations that follow from (l55l are just the 
same as in Minkowski space. Hence we can calculate the Coulomb field as usual. On 
the other hand, the equations of motion receive two changes from the cosmological ex- 
pansion term: the first is that the mass m is now multiplied with the (time-dependent!) 
scale factor a, the second is an additional scalar force induced by a. Note that all 
spacetime dependent functions on the right hand side are to be evaluated at the parti- 
cle's location z{tc), whose fourth component corresponds to etc- Hence, writing out 
all arguments and taking into account that the time coordinate is tc, we have for the 
equation of motion 

z^ = ^jrrr F^a{z)i'' - (-c\'^" + lna,(zO/c) (56a) 

mac[z^ /c) 

= ^07^ - [-cv^^ + i^iVc) <(zO/c)/a,(zO/c) , (56b) 

mac[z^/c) ^ ' 
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where a!^ is the derivative of Oc- 

So far no approximations were made. Now we write = 7(0, u), where v is the 
derivative of z with respect to the conformal time tc, henceforth denoted by a prime, 
and 7 = 1/ y/l — jc? . Then we specialize to slow motions, i.e. neglect effects of 
quadratic or higher powers in vjc (special relativistic effects). For the spatial part of 
(I56bt we get 

z" + z' ia'/ac) = ^(E + z' xB), (57) 
mac 

where we once more recall that the spatial coordinates used here are the comoving 
(i.e. conformal) ones and the electric and magnetic fields are evaluated at the particle's 
position z{tc). 

From the above equation we see that the effect of cosmological expansion in the con- 
formal coordinates shows up in two ways: first in a time dependence of the mass which 
scales with ac, and, second, in the presence a friction term. Let us, for the moment, 
neglect the friction term. In the adiabatic approximation, which is justified if typical 
time scales of the problem at hand are short compared to the world-age (corresponding 
to small 62 in (I36bt ). the time-dependent mass term leads to a time varying radius in 
comoving (or conformal) coordinates of r{tc) oc l/adtc)- Hence the physical radius 
(given by the cosmological geodesically spatial distance), = Ocr, stays constant in 
this approximation. In this way Dicke & Peebles concluded in [49] that electromag- 
netically bound systems do not feel any effect of cosmological expansion. 

Let us now look at the effect of the friction term which the analysis of Dicke & Peebles 
neglects. It corresponds to the decelerating force —va'dac which e.g. for the simple 
power-law expansion d43l becomes 



-V 



t-^ = -vnaf"-^ . (58) 



Clearly it must cause any stationary orbit to decay. For example, as a standard first- 
order perturbation calculation shows, a circular orbit of radius r and angular frequency 
(with respect to conformal time tc) Wc will suffer a relative decay per revolution of 



Ar 

r 



rcvol. Sw, 



Recall that this is an equation in the (fictitious) Minkowski space obtained after rescal- 
ing the physical metric. However, it equates two scale invariant quantities. Indeed, the 
relative length change Ar/r is certainly scale invariant and so is the relative length 
change per revolution. On the right hand side we take the quotient of two quantities 
which scale like an inverse time. In physical spacetime, coordinatized by cosmological 
time t, the right hand side becomes — H/?,uj, where as usual H = a/a and w is the 
angular frequency with respect to t. Hence the relative radial decay in physical space 
is of the order of the ratio between the orbital period and the inverse Hubble constant 
('world-age'), i.e. of order £2 (cf.( l36bt ). 

Since the friction term contributes to the leading-order effect of cosmological expan- 
sion, we conclude that the argument of Dicke & Peebles, which neglects this term, is 
not sufficient to estimate such effects. 
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4.2 Equations of motion in the physical coordinates 



We now show that ^ is indeed arrived at if the friction term is consistently taken into 
account. To see this we merely need to rewrite equation dSTl in terms of the physical 
coordinates given by the cosmological time t and the cosmological geodesic spatial 
distance := a{t)r. We have dtc/dt = 1/a and the spatial geodesic coordinates are 
y := a{t)z. Denoting by an overdot the time derivative with respect to t, the left hand 
side of dSTl becomes 

z" + z' {a'Jac) = a'^--dy. (60) 

This shows that the friction term in the unphysical coordinates becomes, in the physical 
coordinates, the familiar acceleration term ^ due to the Hubble-law. Dividing by a 
equation dSTl and inserting E{z) = Qz/\z\^ and B{z) = 0, we get 



eQ 

m\y\"- 



y-y{d/a) = —j-^y. (61) 



Finally, introducing polar coordinates in the orbital plane we exactly get 



5 Kinematical effects 

It has been suggested in [72] and again in [73] that there may be significant kinemat- 
ical effects that may cause apparent anomalous acceleration of spacecraft orbits in an 
expanding cosmological environment. More precisely it was stated that there is an 
additional acceleration of magnitude He 0.7 • lO^^m/s^, which is comparable to 
the measured anomalous acceleration of the Pioneer spacecrafts. The cause of such an 
effect lies in the way one actually measures spatial distances and determines the clock 
readings they are functions of (a trajectory is a 'distance' for each given 'time'). The 
point is this: equations of motions give us, for example, simultaneous (with respect to 
cosmological time) spatial geodesic distances as functions of cosmological time. This 
is what we implicitly did in the Newtonian analysis. But, in fact, spacecraft ranging 
is done by exchanging electromagnetic signals. The notion of spatial distance as well 
as the notion of simultaneity introduced thereby is not the same. Hence the analytical 
expression of the 'trajectory' so measured will be different. 

To us this seems an important point and the authors of [72] and [73] were well justified 
to draw proper attention to it. However, we will now explain why we do not arrive at 
their conclusion. Again we take care to state all assumptions made. 



5.1 Local Einstein-simultaneity in general spacetimes 

We consider a general Lorentzian manifold {M. , g) as spacetime. Our signature con- 
vention is (+, — , — , — ) and ds is taken to have the unit of length. The differential of 
eigentime is dr = ds/c. In general coordinates {x^}, the metric reads 

ds'^ = g^„dx^ dx" = gttdf + 2gtadt dx"- + gabdx"- dx^ . (62) 
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The observer at fixed spatial coordinates is given by the vector field (normalized to 

g{X,X)=c^) 

^/gttOt 



Consider the light cone with vertex p G A4; one has ds"^ = 0, which allows to solve 
for dt in terms of the dx"- (all functions gab are evaluated at p, unless noted otherwise): 



dti2 = - — dx"± J(^^-^)dx'^dx^ (64) 
9tt V V 9u 9ttJ 

The plus sign corresponds to the future light-cone at p, the negative sign to the past 
light cone. An integral line of X in a neighbourhood of p cuts the light cone in two 
points, (7+ and q^. If tp is the time assigned to p, then tq^ = tp+dti and tq_ = tp+dt2- 
The coordinate-time separation between these two cuts is tq^ — tq_ = dti — dt2, 
corresponding to a proper time ^/gtt{dti — dt2)/c for the observer X. This observer 
will associate a radar-distance dl^ to the event p of c/2 times that proper time interval, 
that is: 

dll = h=( - gab) dx^ dx' . (65) 
V 9tt ) 

The event on the integral line of X that the observer will call Einstein-synchronous 
with p lies in the middle between (7+ and q^. Its time coordinate is in first-order 
approximation given by \{tqj^ + = + \ {dt\ + di2) = + dt, where 

dt := \{dtx + dt2) = -— dx" . (66) 

9tt 

This means the following: the Integral lines of X are parameterized by the spatial 
coordinates {x'^}a=i,2,3- Given a point p, specified by the orbit-coordinates Xp and 
the time-coordinate tp, we consider a neighbouring orbit of X with orbit-coordinates 
Xp + dx"-. The event on the latter which is Einstein synchronous with p has a time 
coordinate tp + dt, where dt is given by (I66t . or equivalently 

e ■=dt + ^dx'' = 0. (67) 
9tt 

Using a differential geometric language we may say that Einstein simultaneity defines 
a distribution 9 = 0. 

The metric d62b can be written in terms of the radar-distance metric h d65b and the 
simultaneity 1-form 6 as follows: 

ds^ = g^^.dx'' dx" =gud^-h, (68) 

showing that the radar-distance is just the same as the Einstein-simultaneous distance. 
A curve 7 in 7W intersects the flow lines of X perpendicularly iff 6* (7) = 0, which is 
just the condition that neighbouring clocks along 7 are Einstein synchronized. 
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5.2 Application to isotropic cosmological metrics 



We consider isotropic cosmological metrics. In what follows we drop for simplicity 
the angular dimensions. Hence we consider metrics of the form 

ds'^ = c^dt^ - a{tfdr^ . (69) 

The expanding observer field is 

ot 

The Lagrangian for radial geodesic motion is L = ^ (c^P — a^r^) , leading to the only 
non-vanishing Christoffel symbols 

^rr ~ ~9" ' ^tr = — ='■ H . (71) 

a 

Hence X is geodesic, since 

VxX = T^tdf, = 0. (72) 

On a hypersurface of constant t the radial geodesic distance is given by ra{t). Making 
this distance into a spatial coordinate, r*, we consider the coordinate transformation 

1 1— > := t , r 1-^ := a{t)r . (73) 

The field d/dt^ is given by 

^ 9 TT 9 
= Hr — . (74) 

c?t* dt dr 

In contrast to dTOb . whose flow connects co-moving points of constant coordinate r, the 
flow of (I74t connects points of constant geodesic distances, as measured in the surfaces 
of constant cosmological time. This could be called cosmologically instantaneous 
geodesic distance. It is now very important to realize that this notion of distance is not 
the same as the radar distance that one determines by exchanging light signals in the 
usual (Einsteinian) way. Let us explain this in detail: 

From (1731 we have adr = dr^ — r^Hdt, where H := a/a (Hubble parameter). Rewrit- 
ing the metric ( l69l in terms of and r* yields 

= 0^(1 - {Hrjcf) dtl - drl + 2Hr^ dt dr^ 
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Hence the differentials of radar-distance and time-lapse for Einstein-simultaneity are 
given by 
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Let the distinguished observer (us on earth) now move along the geodesic = 0. 
Integration of (1761 from = to some value then gives the radar distance as 
well as the time lapse At* as functions of the cosmologically simultaneous geodesic 
distance r*: 

h = {c/H) sm-\Hn/c) ^ n{l + ^{Hr,/cf + 0{3)} (77a) 
At, = {1/2H) ln(l - (Hn/cf) ^ {n/c){-^{Hn/c) + 0{2)} (77b) 

Combining both equations in dTTt allows to express the time-lapse in terms of the 
radar-distance: 

At, = H-^ ln(cos(F/,/c)) ^ {h/c){-liHk/c)+0{2)} . (78) 

Now, suppose a satellite S moves on a worldline (t*) in the neighbourhood of our 
worldline = 0. Assume that we measure the distance to the satellite by radar 
coordinates. Then instead of the value we would use /* and instead of the argument 
t* we would assign the time t* — At* which corresponds to the value of cosmological 
time at that event on our worldline that is Einstein synchronous to the event (t*, r*); 
see Figure|2l 

Hence we have 

= (c/H) sm~^{n{U + AU)H/c} (79a) 
^ r,-^{v/c){Hc){r,/cf , (79b) 

where (I79bt is (I79al ) to leading order and all quantities are evaluated at t*. We set 

V = r*. 

To see what this entails we Taylor expand in t* around t* 
choice): 

r*(t*) = ro + uot* + |aot* H 

and insert in ( I79bt . This leads to 

/*(t*) = ro + woi* + ^aotl H , (81) 

where, 

fo = ro-(//c) i(^;o/c)(ro/c)2 (82) 

vo = vo- (He) {vo/c)\ro/c) (83) 

ao = ao - (He) {{vo/cf + (ro/c)(^;o/c)(ao/c)} (84) 

These are, in quadratic approximation, the sought-after relations between the quantities 
measured via radar tracking (tilded) and the quantities which arise in the (improved) 
Newtonian equations of motion (not tilded). 

The last equation (l84l i shows that there is an apparent inward pointing acceleration, 
given by He times the {v/c)^ + ■ ■ ■ term in curly brackets. The value of He is indeed 
of the same order of magnitude as the anomalous Pioneer acceleration, as emphasized 
in [72, 73]. However, in contrast to these authors, we do get the additional term in 



(just a convenient 
(80) 
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us (r» = 0) 




satellite (r*(i*)) 



•■■■'if* = const. 



Figure 2: The observer ('us') moves on the geodesic worldline r^, = and measures 
the spatial distance to a satellite by exchanging electromagnetic signals (radar coor- 
dinates). The satellite's worldline is analytically described by a function r*(t^,). The 
surface of constant cosmological time (tilted dashed line) intersects our worldline 
at A and the satellites worldline at B. The coordinate corresponds to the geodesic 
distance in that hypersurface of constant cosmological time, i.e. to the length of AB as 
measured in the spacetime geometry. However, using radar coordinates, the observer 
defines B to be simultaneous to his event A' and attributes to it the distance l^. Hence 
instead of r^,(t*) he uses ^*(r^,(t* — At^)), where l^:{r^) is given by (I77at . This leads 
todllli. 
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curly brackets, which in case of the Pioneer spacecraft suppresses the He term by 
13 orders of magnitude! Hence, according to our analysis, and in contrast to what is 
stated in [72, 73], there is no significant kinematical effect resulting from the distinct 
simultaneity structures inherent in radar and cosmological coordinates. We should 
stress, however, that this verdict is strictly limited to our interpretation of what the 
kinematical effect actually consists in, which is most concisely expressed in (|7917 

6 Summary and outlook 

We think it is fair to say that there are no theoretical hints that point towards a dynam- 
ical influence of cosmological expansion comparable in size to that of the anomalous 
acceleration of the Pioneer spacecrafts. There seems to be no controversy over this 
point, though for completeness it should be mentioned that according to a recent sug- 
gestion [70] it might become relevant for future missions like LATOR. This suggestion 
is based on the model of Gautreau [14] which, as already mentioned in Section |3j we 
find hard to relate to the problem discussed here. Rather, as the (a/a)-improved New- 
tonian analysis in Section [Altogether with its justification given in Sections |3] and |3 
strongly suggests, there is no genuine relativistic effect coming from cosmological 
expansion at the levels of precision envisaged here. 

On the other hand, as regards kinematical effects, the situation is less unanimous. It is 
very important to unambiguously understand what is meant by 'mapping out a trajec- 
tory', i.e. how to assign 'times' and 'distances'. Eventually we compare a functional 
relation between 'distance' and 'time' with observed data. That relation is obtained 
by solving some equations of motion and it has to be carefully checked whether the 
methods by which the tracking data are obtained match the interpretation of the coor- 
dinates in which the analytical problem is solved. In our way of speaking dynamical 
effects really influence the worldline of the object in question, whereas kinematical ef- 
fects change the way in which one and the same worldline is mapped out from another 
worldline representing the observer. 

The latter problem especially presents itself in a time dependent geometry of space- 
time. Mapping out a trajectory then becomes dependent on ones definition of 'simul- 
taneity' and 'simultaneous spatial distance' which cease to be unique. An intriguing 
suggestion has been made [72, 73] that the PA is merely a result of such an ambiguity. 
However, our analysis suggests that no significant relativistic effects result within the 
Solar System, over and above those already taken into account, as e.g. the Shapiro 
time delay. 

What has been said so far supports the view that there is no interesting impact of 
cosmological expansion on the specific problem of satellite navigation in the Solar 
System. However, turning now to a more general perspective, the problem of how 
local inhomogeneities on a larger scale affect, and are affected by, cosmological ex- 

' Our equation <78t corresponds to equation (10) of [72]. From it the authors of [72] and [73] im- 
mediately jump to the conclusion that there is "an effective residual acceleration directed toward the 
centre of coordinates; its constant value is He". We were unable to understand how this conclusion is 
reached. Our interpretation of the meaning of <79> does not support this conclusion. 
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pansion is of utmost importance. Many scientific predictions concerning cosmological 
data rely on computations witiiin the framework of the standard homogeneous and 
isotropic models, without properly estimating the possible effects of local inhomo- 
geneities. Such an estimation would ideally be based on an exact inhomogeneous 
solution to Einstein's equations, or at least a fully controlled approximation to such 
a solution. The dynamical and kinematical impact of local inhomogeneities might 
essentially influence our interpretation of cosmological observations. As an example 
we mention recent serious efforts to interpret the same data that are usually taken to 
prove the existence of a positive cosmological constant A in a context with realistic 
inhomogeneities where A = 0; see [47] and [58]. 

To indicate possible directions of research, we stress again that there are sev- 
eral approaches to the problem of how to rigorously combine an idealized local 
inhomogeneity — a single star in the most simple case — with an homogeneous and 
isotropic cosmological background. We mentioned that of Einstein & Straus [7] and 
its refinement by Schiicking [34], that of Gautreau [14], that by Bonnor [44, 45, 46, 3], 
and especially the classic work by McVittie [59] that was later elaborated on by 
Hogan [19] and properly interpreted by the penetrating analysis of Nolan [30, 31, 32]. 
Nolan also showed that only in the case k = does McVittie's solution represent a 
central mass embedded into a cosmological background. The problem is that, due to 
the non-linearity of Einstein's equations, the spherical inhomogeneity does not show 
just as an addition to the background. Hence a notion of quasi-local mass has to be 
employed in order to theoretically detect local mass abundance. However, as is well 
known, it is a notoriously difficult problem in General Relativity to define a physi- 
cally appropriate notion of quasi-local mass. Workable definitions only exist in special 
circumstances, as for example in case of spherical symmetry, where the concept of 
Misner-Sharp mass can be employed, as explained in Section l3!n 

As a project for future research we therefore suggest to further probe and develop 
applications of the spatially flat (k = 0) McVittie solution, taking due account of 
recent progress in our theoretical understanding of it. As a parallel development, the 
implications of Gautreau's model should be developed to an extent that allows their 
comparison with those of McVittie's model. 
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A Additional material 



In this section we collect some background information which was implicitly used 
throughout the text. 



A.l A simple estimate of the dynamical effect of cosmological expansion 

The radial acceleration due to the Newtonian Sun attraction is given by: 

The geodesic distance, r*, between two freely falling bodies in an expanding FRW 
universe d37b . measured on a hypersurface of constant cosmological time, varies in 
time according to the Hubble-law = i7r*. The related 'acceleration' is then f^, = 
fifH+r^H = r^{H'^+H) = r^a/a = —qH^r^. At the present time (see Section lA.3ll 
we have (suppressing the asterisk): 

r |cos„..acc. = - r pa 4 (r/10 AU) W~^^ m/s" . (86) 
a 

This naive derivation of the dynamical effect of the cosmological expansion may be 
confirmed by the fully general relativistic treatment, as showed in Section |3l Notice 
that, according to our measurements, the universe is presently in a phase of accelerated 
expansion, hence (l86b results in an acceleration pointing away from the Sun. This is 
in the opposite direction of the Pioneer-effect d92t and also smaller by 14 orders of 
magnitude. 



A.2 The Schiicking radius for different astronomical scales 

In the following we evaluate the radius of the Schiicking vacuole (flOl for various 
characteristic central masses. Using as cosmological matter density pm = QmPc (see 
Appendix IA.3t . we get: 

• Solar System scale: rs{MQ) = 570 ly, which is much larger than the average 
distance between stars in the Milky Way, being about 10 ly. 

• Galaxy scale: rsiMMw) = 3 — 4 Mly, which is again too big since this would 
also include other galaxies such as the Large and the Small Magellanic Cloud, 
as well as several dwarf galaxies. 

• Cluster scale: ^^(Mlg) = 5 — 7 Mly, which is just about the threshold since 
the nearest object not belonging to this cluster is NGC 55 (a galaxy belonging 
to the Sculptor Group) about 5 Mly away. 

• Supercluster scale: r5(Mvsc) = 57 Mly, which is inside the Virgo Superclus- 
ter whose radius is about 100 Mly. 

This shows that the Einstein-Straus matching works at best from and above cluster 
scale. 
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A.3 Some astronomical and cosmological data 

For the convenience of the readers, we collect some relevant numerical information. 
Length units 

1 AU = 149.6 • 10<^ km = 1.5 • 10^^ m = 492 Is = 8.2 hnin = 1.58 • 10"^ ly 
1 pc = 3.26 ly 

Time units 

lyr = 3.1-10^s 

Age of the Universe ^ 13.7 • 10^ yr = 4.32 • 10^^ s 

Velocity units 

1 AU/yr = 4.74 km/s 

Mass units 

1 Mq = 2 • 10^° kg (= sometimes referred to as twice the mass of the average star in 
the Milky Way) 

The Universe 

Our galaxy: Milky Way (MW) 
Number of stars in the MW: 2-4 • 10^^ 
Mmw = 2-4 • IO^^Mq 
Diameter: 100 kly 

Average distance between stars in the MW: 10 ly 

Nearest galaxy: Large Magellanic Cloud. Mass: 10^^ Mq, distance: 170 kly, diame- 
ter: 30 kly 

Our cluster: Local Group (LG) 

Number of stars in the LG: 7 • 10^-*^ 
Mlg = 7-20- 1O"M0 
Diameter: lOMly 

Nearest clusters: Sculptor Group (distance: 10 Mly) and Maffei 1 Groups (distance: 

lOMly). 

Nearest galaxy: NGC55 (Sculptor Group), distance: 5 Mly. 

Our supercluster: Virgo Supercluster (VSC) 
Number of stars in the VSC: 2 • 10^^ 

Mvsc = 1 • W^^Mq 
Diameter: 200 Mly 

Nearest supercluster: Centaurus (distance: 200 Mly) 

Nearest cluster: A3526 (Centaurus Supercluster), distance: 142 Mly 
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Cosmology data 

Hubble parameter (today): Hq = ho ■ 100 km s"^ Mpc"^ = /io/(3.08 • lO^^s), 
where ho = 0.7 

Critical density (today): = 3H^/8ttG = hi ■ 1.89 • lO'^^ g/cm^ 
Definition of the cosmological parameters: 

^m'-= Pml Pc, (87) 

•= Pk/Pc = , (88) 

■■= -kc'/iHoaof . (89) 

One has Qjn + + = ^ (cosmological triangle), where todays values are given 
by: 

(17™, J^A, Qfc) (1/3,2/3,0). (90) 
Thus from qo = (1/2)17^ — it follows that qq —1/2. 



Pioneer 10 and 11 data 



The following data are taken from [75]. 





PIO 


Pll 


Launch 


2 Mar 1972 


5 Apr 1973 


Planetary 


Jupiter: 4 Dec 1973 


Jupiter: 2 Dec 1974 


encounters 




Satum: 1 Sep 1979 


Tracking data 


3 Jan 1987 - 22 Jul 1998 


5 Jan 1987 - 1 Oct 1990 


Distance from the Sun 


40 AU - 70 AU 


22.4 AU- 31.7 AU 


Light round-trip time 


11 h- 19h 


6h-9h 


Radial velocity 


13.1 Km/s- 12.6 Km/s 


N.A. 



Table 1: Some orbital data of the Pioneer 10 and 1 1 spacecrafts. 



Tracking system 

Uplink frequency as received from Pioneer (approx.): Uu^2 = 2.11 GHz. 
Downlink frequency emitted from Pioneer: f,^ 2 = 7^J^u,2 = 2.292 GHz. 
Spacecraft transponder turnaround ratio: T = 240/221. 

Measured effect 

Measured is an almost constant residual (meaning after subtraction of all the known 
effects) frequency drift of the received tracking signal. The drift is a blue-shift at the 
constant rate 

2/ = (5.99 ±0.01)10-'' Hz/s (91) 

which, if interpreted as a special-relativistic Doppler shift, can be rewritten as an ac- 
celeration pointing towards the Earth (or Sun) of modulus 

ap = (8.74 ± 1.33) 10"^° m/s^ . (92) 
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